iju- ir-iiui 



DESY 06-167 



On the equivalence of Reggeon field theory in zero transverse 
dimensions and reaction-diffusion processes 



S. Bondarenko a ) * , L. Motyka^ f , A.H. Mueller d > * , A.I. Shoshi e ) § and B.-W. Xiao d ) 1 



a ) // Institute for Theoretical Physics, University of Hamburg, 22761 Hamburg, Germany 
b ) DESY Theory Group, 22603 Hamburg, Germany 
c ) Institute of Physics, Jagellonian University, 30-059 Krakow, Poland 

d ) Physics Department, Columbia University, New York, NY-10027, USA 

e ) Fakultdt fur Physik, Universitat Bielefeld, D-33501 Bielefeld, Germany 

September 20, 2006 



Abstract 

The Reggeon field theory in zero transverse dimensions is investigated. Two versions of the 
theory are considered: one that allows at most triple pomeron interactions and the other that 
embodies an additional 2 — > 2 quartic Reggeon coupling. The behavior of the scattering amplitude 
at asymptotic rapidities is obtained in both cases. In an s-channel picture of the high energy 
scattering both models can be viewed as reaction-diffusion processes. We derive known results 
in Reggeon held theory rather easily using the reaction-diffusion formalism. We find that some 
results which are surprising from the Reggeon field theory point of view turn out to have a simple 
interpretation from the reaction-diffusion point of view. 
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1 Introduction 



The high energy scattering in QCD may be described in terms of a non-local field theory of QCD 
pomerons. In this formulation the basic degrees of freedom — the BFKL pomerons ^[21131 13 — merge 
and split according to non-local multi-pomeron vertices [SJ H3 [7| |B1 H3 EEB EEU dH E3 • This approxima- 
tion is consistent with the unitarity of the S'-matrix. Alternative approaches to the unitarization of 
high energy amplitudes use extensively methods of statistical physics combined with field-theoretical 
calculations in QCD, see e.g. Refs. |14|-|34|. The coordinate space in all those formulations is the 
transverse plane to the axis of the high energy collision. In all cases the emerging evolution equa- 
tions for high energy amplitudes are non-local functional equations which are, so far, prohibitively 
difficult to solve exactly. Therefore, it is useful to study a substantially simplified model of the inter- 
acting pomeron theory, in which the dependence of the pomeron fields on transverse coordinates is 
discarded. The emerging simplified model - Reggeon field theory (RFT) in zero transverse dimensions 
- was proposed and formally solved long time ago |351 1361 1371 1381 1391 14U| . This toy model exhibits 
some features which are present also in the more realistic QCD pomeron field theory. Therefore, RFT 
in zero transverse dimensions offers a useful testing ground for computational methods and recently 
it enjoys revived interest |411 1421 1431 1441 1451 146j . In this paper we will explore the equivalence be- 
tween RFT and the statistical reaction-diffusion picture of high energy scattering that holds in zero 
transverse dimensions. 

In Section 2 we consider the zero dimensional Reggeon field theory. Equations of RFT are derived 
from the Langevin formulation and from the Lagrangian of RFT. Two distinct realizations of RFT are 
studied: a minimal RFT with triple pomeron vertices only (which will be abbreviated as MRFT) and a 
theory with an additional 2 — ► 2 pomeron interaction which is equivalent to a simple reaction-diffusion 
s-channel model (which will be called the reaction-diffusion RFT or RD-RFT). Asymptotic behavior 
of scattering amplitudes in both RFTs is obtained. The minimal RFT yields a counter-intuitive 
exponential decrease for the scattering amplitude T(Y) ~ exp(— EqY) at very large rapidities Y, 
with the leading order result of Eq given by Eq ~ exp(— /i 2 /2A 2 ) where /x denotes the pomeron 
intercept and A the triple pomeron coupling. Clearly, the dependence of Eq on the triple pomeron 
coupling is non-perturbative. 

The presence of a quartic 2 — ► 2 pomeron coupling, A', changes this picture. Namely, since Eq = 
for A' = A 2 ///, the scattering amplitude approaches a constant value at asymptotic rapidities instead 
of going to zero. In this case, the exact asymptotic solution of the RFT is given. Furthermore, 
we briefly discuss two possible options of multi-pomeron couplings to elementary external sources. 
Firstly, the eikonal couplings that were usually assumed in RFT. In the other scenario only a single 
pomeron coupling to the elementary external source is allowed. The latter scheme is realized in QCD 
where (at the leading logarithmic approximation) only single BFKL pomeron couples to an elementary 
colour dipole [HI El El El • It is also favored by the s-channel reaction-diffusion picture in which each 



2 



elementary particle in the cascade may couple to at most one f-channel exchange. In this context, it 
is particularly interesting to consider the RFT counterpart of the reaction-diffusion model in which 
only single coupling to external sources are allowed with the strength a <C 1 related to the pomeron 
intercept \x = a and the triple pomeron coupling A = a 2 . In this model the high energy asymptotics, 
Y — > oo, of the scattering amplitude happens to be close to the unitarity limit: 

which should be contrasted with the result of RD-RFT with the same parameters and eikonal couplings: 

T (Y) - f^k- (1.2) 

In Section 3 we give a complementary analysis of the scattering amplitudes in zero transverse 
dimensions in the probabilistic reaction-diffusion framework. It turns out, that the minimal RFT may 
be represented in the s-channel picture, analogously to the reaction-diffusion process. The transition 
rates, however, do not obey the constraints of the Markovian process - the transition coefficients 
are negative that prohibits the genuine probabilistic interpretation. Nevertheless, the quasi-reaction- 
diffusion picture provides an efficient calculational framework. Thus, we obtain results for high rapidity 
asymptotics of scattering amplitudes that coincide with ones from the conventional RFT framework, 
but the interpretation is more transparent. We also demonstrate the equivalence of the results obtained 
within the proper reaction-diffusion framework with the results of the corresponding RFT. Finally, 
the representation of the eikonal model in the reaction-diffusion framework is given. Curiously, the 
"grey disc" limit of the scattering amplitudes given by Eq. (|l.'2jl emerges as a consequence of the 
normalisation of the scattering states that is smaller than one. 

Finally, we note that a different s-channel picture has recently been elaborated in Ref. |46| which 
is very different than the one we have here discussed. The eikonal couplings are used and boost 
invariance is preserved by having arbitrarily high Reggeon vertices in the hierarchy equations. The 
system does not reach a fixed point, but rather reaches a steady rate of increase of particles with 
increasing rapidity. It will be interesting to try and further clarify which of these pictures, if either, 
has the closest relationship with four-dimensional QCD. 



2 Reggeon field theory in zero transverse dimensions 

2.1 Reggeon field theory from the Langevin equation. 

The general form of the Langevin equation in the Ito formulation is 

du{y + dy) = D^uiy^dy + D 2 (u{y)) duj{y) , (2.1) 
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where u(y) is a stochastic function of an evolution variable y, dy is an infinitesimally small step of y, 
and to{y) is related to a random noise variable v(y), 

My) = " (j/) d y ■ (2-2) 

The noise is taken to be Gaussian, 

(i/(y)> = o, MvMi/)> = (2-3) 

Let us consider a dynamical quantity G(it) with it distributed with the probability density p(y, u). 
It follows from the Ito calculus that the average value of this variable 



(G) y = J dup(y,u)G(u), (2.4) 
evolves according to the following master equation: 

_/ DlW ^\ + (DMim\ . (2 . 5) 



dy \ du / y \ 2 du 2 

Thus, assuming that moments of the variable u 



(u k ) y = / d^(y, U ) U fc (2.6) 



obey the hierarchy of (deterministic) equations 



dy 



k ( Lq^K ) + - 1) ( ^^u k - z ) . (2.7) 



which closes if D\(u) and D^u) are analytic functions of u. With an appropriate choice of the 
functions D\(u) and D%(u) this hierarchy of equations may be made equivalent to the equations of 
the Reggeon field theory in zero transverse dimensions. We choose the following form of D\[u) and 
D 2 (u) 

Di(u) = au-pu 2 , (2.8) 
£>!(«) = 2 (a/ if - /?'n 2 ), (2.9) 

which is equivalent to RFT with triple pomeron vertices and a quartic 2^2 pomeron vertex. Indeed, 
with (|2,8|) and (|2.9|) . one obtains from ()2.7|) a closed set of equations 

d{ - U p y = ak{u k ) y + a'k(k-l){u k - 1 ) y - (3k{u k+l ) y - f3' k{k— 1) (u k ) y . (2.10) 
ay 

As it will be shown, those equations are the same as equations of RFT if the /c-pomeron amplitude 
3>fc(y) is connected to the fc-th moment (u k ) y in the following way, 

*„(3,) = -(-«) fc ^, (2-H) 
and the parameters a' , (3 and 0' are suitably adjusted. 
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2.2 Hamiltonian formulation of RFT 



The analysis in this and the following parts of this section recapitulates the formulation and the key 
results that were obtained within RFT long ago [SHI EH EH EM EH • We describe them in some details 
in order to facilitate a comparison with the following analysis within the s-channel framework. 

Originally, RFT was formulated as a field theory (or as quantum mechanics in zero transverse 
dimensions) of pomerons. The basic degrees of freedom in this formulation are the Gribov fields ij) 
and ip + that create and annihilate the pomeron. The action defining the theory with triple pomeron 
couplings only (MRFT) is defined in the following way: 

S= [ dy {ip + d y ^j - fi%b + %b + (V> + + ip)ip }, (2.12) 



where \x is the bare intercept of the pomeron and A is the coupling of the triple pomeron interaction. 
Note that the action is symmetric under the transform 

^ <— » V + , y -» -v, (2-13) 

which corresponds to the symmetry between the target and the projectile. After redefinition of the 
Gribov fields, q = iip + and p = tip the action may be rewritten in terms of real quantities 

S = J dy {pd y q + fiqp - \q(q + p)p } . ( 2 -14) 

Then, the Hamiltonian of the problem is given by 

H = fiqp — Xq(q + p)p . (2-15) 

The commutation relation of the pomeron annihilation and creation operators, [/!/>, = 1, implies 
[PjQ] = — 1 which may be realized for instance by identification of q with a position operator and p 
with a differential operator, 

P = - 7r . 2.16 
oq 



In this representation the Hamiltonian takes the form 

H = -{^q- \q 2 )^- -Xq^-r. (2.17) 



— - A — 

dq ^ dq 2 



In the following the ratio g = ji / A will be used, which is the large parameter of our model. 

The interacting pomeron system whose state is given by *&(y, q) evolves according to a Schrodinger 
equation (with an imaginary evolution variable) 

-^- = iTtf. (2.18) 
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The complete definition of the problem requires that the domain of q and the boundary conditions 
are specified. This, however, depends on the particular choice of the couplings. For the minimal RFT 
it was shown that it is sufficient to consider < q < oo with ^f(y, q) that obeys the conditions 

= 0, Km d q V(y,q) = (2.19) 

at all rapidities y. 

The state ^(y,q) embodies the complete information about /c-pomeron amplitudes &k(y) in the 
i-channel, 

oo 

*(y,q) = £ $k(y)q k - (2.20) 

fe=l 

The following hierarchy of equations, equivalent to (|2.18l) . 

— * = nk® k + \k{k + l)$ k+1 - A(fc-l)* fc -i, k = l,2,... (2.21) 
ay 

is identical to the hierarchy of equations (|2.1()j) describing the evolution of moments (u k ) in the 
Langevin formulation after the substitution (|2.1H) provided that it is set 

a = [i, aol = p/a = A, p' = 0. (2.22) 

Note, that the target-projectile symmetry of RFT imposes the constraint 

P = a 2 a' (2.23) 

on the Langevin formulation. 



2.3 External couplings and scattering amplitudes 

The scattering amplitudes in the Reggeon field theory may be obtained after the bulk action given by 
Eq. (|2.12l) is supplemented by the action describing the coupling of the pomeron system to external 
particles. In the original formulation of RFT the eikonal couplings were assumed, that is the coupling 
of k pomerons to the external source was given by a product of k single pomeron couplings multiplied 
by the combinatoric symmetry factor, l/k\ for k outgoing pomerons 1 . In this scenario, the initial 
condition of the RFT Hamiltonian evolution is obtained by a resummation of multi-pomeron source 
amplitudes: 

oo <x> ( \k k 

HV = 0, g) = £ $k(0)q k = £ ~ [ ~ 9 ,y Q = 1 - exp(- 5l g), (2.24) 

k=i fc=l fc - 

1 This assumption seems to be rather natural for interacting point-like objects, but it may be not applicable in the case 
of QCD pomerons. In particular, due to a complex non-local nature of the BFKL pomeron, only one BFKL pomeron 
can couple to an elementary colour dipole in the leading logarithmic approximation. 
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where g\ is the coupling of the pomeron to the source. Thus, the scattering amplitude in RFT at 
rapidity y is given by 

oo 

T(Y; 9l ,g 2 ) = Y,$k(Y)gl = *(Y,g 2 ), (2.25) 
fc=i 

where g 2 is the pomeron coupling to the projectile. 

In the alternative realization of RFT only a single pomeron coupling is possible to an elementary 
(thus not evolved) external source. This assumption naturally emerges if one considers dipole-dipole 
scattering in QCD. Then, the initial condition of the Hamiltonian evolution is given by: 

*(y = 0, g ) = *i(0)g = giq, (2.26) 

and the scattering amplitude may be obtained as 

d*{Y,q) 



T(Y; gi ,g 2 ) = g 2 $ 1 (Y) = g 2 



dq 

2.4 High rapidity asymptotics of the minimal RFT 



(2.27) 

<J=0 



The solution to the Hamiltonian problem (|2.18|) may be given in terms of an infinite summation over 
a discrete set of eigenstates with wave function tp n (q) 

oo 

*(y,q) = ]T \ n e- E ^Ml) , (2-28) 

n=0 

with the eigenvalues of the Hamiltonian H denoted by E n , n = 0,1,2,.... Clearly, the asymptotic 
rapidity dependence of the scattering amplitude is governed by the lowest eigenvalue Eq. 

Minimal RFT is known to exhibit a counter-intuitive feature of an asymptotic decrease of the 
scattering amplitude down to zero at Y — > oo. This happens because the lowest eigenvalue Eq of 
the RFT Hamiltonian is positive. We shall recall a heuristic derivation of Eq, given in Ref. [HHl that 
demonstrates this explicitely. 

First, let us transform the eigenfunction corresponding to Eq: 

Mq) = e^ 2 / 4 f (q) , (2.29) 
that leads to the following equation for fo(q) : 

H'fo(q) = f-J^ " \ + \{q ~ Q?^ fo(q) = f|/o(9), (2-30) 
and the boundary conditions in Eq. (|2.19jl impose the following behavior of fo(q): 

fo(q = 0) = (2.31) 
f (q -» cxd) cx e -? 2 /4 + W /2 _ (2 32) 
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This eigenvalue problem may be solved approximately starting from a simpler problem of the ground 
state of the harmonic oscillator, with only modification due the initial condition Eq. (|2.31[) . In this 
case, the Coulomb term Eo/Xq is treated as a small perturbation. Then the value Eq must be adjusted 
so that the lowest eigenvalue of the operator H' — Eo/Xq vanishes. Fortunately, the eigenfunctions 
of operator H' in Eq. (EOnjl defined on L|(0, oo) are known. These are parabolic cylinder functions 
D v (q - g): 

H'D u {q-g) =uD„(q- e ). (2.33) 

It emerges self-consistently that v quickly vanishes with increasing g, specifically v ~ e~ g2 / 2 . The 
initial condition at q — > implies 

D v {-g) = 0, (2.34) 

fulfilled by v = vq that may be obtained approximately after performing asymptotic expansion of 
D u (—g) in terms of the large parameter g: 

vq ~ e e -e 2 l\l + 0(l/g 2 )). (2.35) 



In the first order of the perturbative analysis in Eq, the contribution of the Coulomb term is given 

This correction must be adjusted so that it cancels to zero the energy of the operator H' — E^/Xq . 
Thus it follows from Eq. (t2~33l and Eq. dPfll that 

Eq 



uq, (2.37) 



leading to the energy of the ground state 



E ~ -^fe-^/^l + Oil/g 2 )), (2.38) 



and the behavior of the RFT wave function at large rapidities, 

*(y,q) ~ ^ 2 / i D I/0 (q-g)e-^^ ee ~ e2/2 . (2.39) 

Therefore, the scattering amplitudes at fixed external couplings decrease exponentially with increasing 
rapidity 



T(Y,g 1 ,g 2 ) ~ exp 



Y „2 



(2.40) 



f 2ir 

2 

Note, that the exponent Eq is exponentially small at large g, Eq ~ e~ e , and that the splitting of 
Eq from zero is a non-perturbative effect in the triple pomeron coupling, A = fi/g. This exponential 
decrease of the amplitude may be interpreted as a result of a tunneling phenomenon |36 | I39 [ HO] . In 
fact, the tunneling rate was explicitly calculated using the instanton calculus |4Uj which gives exactly 
the y-dependence given by Eq. (|2.40|) . Thus, the semi-classical calculation of the (eikonal) MRFT 
scattering amplitude at large Y yields: 



T(Y, gi ,g 2 ) ~ (1 - e^i) (1 _ e -ew) eX p 



Y „2 



lige-fl 2 



(2.41) 
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2.5 Markovian RFT 



Another special case of RFT is obtained for the choice 

a' = a, = (3 = a 3 , (2.42) 

in the Langevin equation. In this case the evolution of the system may be represented in terms of 
a Markovian reaction-diffusion process in O-dimensions with the 1 — ► 2 particle splitting probability 
a and 2 — > 1 particle merging probability /3 = a 3 . To be precise, in this case hierarchy of evolution 
equations 1)2. 10|) for moments (u k ), 



d{u k 



^ = ak(u k ) y + ak(k-l)( u^ 1 ) y - a 3 k( u k+1 ) y - a 3 k (k - 1) ( u k ) y , (2.43) 



dy 

is identical to the hierarchy of evolution equations (|3.11|) for factorial moments n,( k ' of particle number 
in the reaction-diffusion process. This version of RFT, which we denote as RD-RFT, is defined by the 
following hierarchy of evolution equations for /c-pomeron amplitudes 

^ = fik^ k + \k(k + l)<S> k+1 - A(fc-l)* fc _i - \'k(k + l)<S> k (2.44) 

with the pomeron intercept // = a, the triple pomeron coupling A = a 2 and the quartic 2^2 pomeron 
coupling A' = a 3 . Consequently, the Hamiltonian of the system reads 



f) f) 2 



H M = {\q 2 - m )— + (\'q 2 - \q)—. (2.45) 



In general, the Hamiltonian H m factorizes into a g-dependent part (without derivatives) and a d/dq- 
dependent part provided that A'//i = (X/fi) 2 . This condition is obeyed in particular by values of A 
and A' imposed by the relation to the reaction-diffusion model and the Hamiltonian reads 



H 



M 



dq dq 2 



In order to find the scattering amplitude we consider the quantum mechanical problem defined by this 
Hamiltonian: 

- — = H M V . (2.47) 
dy 

The boundary conditions that should be imposed on *S>(y,q) in the regular singular points q = and 
q = g = 1/a are the following: 

*(y, g = 0)=0 (2.48) 
fy(y, q) is an analytic function around q = 1/a. (2.49) 

The solution may be represented in terms of the eigenfunctions i/j n (Q) and eigenvalues of the Hamil- 
tonian Hj[j 

n=oo 

*(y,q) = Yl Ke- E " y Mq)- (2-50) 

n=0 
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Clearly, the state characterized by 

E = (2.51) 

Ml) = 1 " e~ ?/a , (2.52) 

is an eigenstate. In fact, this is the ground state of Hm which dominates the scattering amplitude 
at high rapidities y. All the higher states with n = 1,2, .. . obey the condition 2 ip n (g) = 0, see e.g. 
Ref. 03]. Therefore, the weight Ao in Eq. (|2.50f) may be found from the initial condition at y = 0, 

q,( y = 0,q = 1/a) = A Vo(l/a) = A (1 - e" 1 ^ 2 ), (2.53) 

and the solution at large y is approximately given by 

* ^i% ^-e- q/a ). (2.54) 

Assuming the eikonal couplings to the external sources, characterized by coupling constants g\ 
and <72> the initial condition takes the form ^(0, q) = 1 — exp(—g\q), and the scattering amplitude at 
high Y reads 

(I - e -9i/a) (l _ e -92/a\ 

T(Y, 9l ,g 2 ) ~ U \_ } e [ _ 1/a2 >-. (2.55) 

For g\ = §2 = a -C 1 and Y — > oo this amplitude tends to (1 — l/e) 2 /(l — e" 1 /" 2 ) which is far from 
the unitarity limit. 

In the standard reaction-diffusion picture only coupling of single pomeron to the elementary ex- 
ternal source is allowed, with the strength a, and g\ = gi = a. Then, with the initial condition 
9(0, q) = q/a, one obtains for aY ^> 1 



f(Y, gi ,g 2 ) 



(gi/a) (52 /a) 



1 — p~ x l a 



! (2.56) 



91=92=0 1 — e 1 / q2 

As we will see in the next section, the results for the amplitude given by Eq. (|2.55|) or Eq. (|2.56j) 
may be reproduced in the reaction-diffusion model confirming the equivalence of the approaches. 



3 High energy behavior of the scattering amplitude in diffusion 
reaction framework 

In this section we use the reaction-diffusion formalism to reproduce the above results in Reggeon field 

theory. This yields two benefits. First the results of Reggeon field theory emerge rather easily using 

the reaction-diffusion formalism and, secondly, some seemingly surprising results from the Reggeon 

field theory point of view have a simple interpretation from the reaction-diffusion point of view. 

2 This may be easily seen from an analysis of the power-series solutions of the eigenvalue equation, HmiPu = E n tp n , 
around the regular singular point q = g. Namely, for any E n 7^ there exists an analytic solution which approaches 
zero as q — g at q — » g. The other independent solution does not vanish at q = g but it contains non-analytic terms 
~ log(f> — q) and thus it is rejected. 
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3.1 3P coupling only 



High energy scattering where only a three pomeron interaction is present in the Reggeon field theory 
is a very natural theory to study from the t-channel (Reggeon field theory) perspective. It does 
not naturally correspond to a reaction-diffusion theory but, formally, one can force it into such a 
description. As we have seen in the previous section the hierarchy equation for the factorial moments, 



n v * is 



dy 
with 



n<*> 



= akn^ + ak(k - l)^*" 1 ) - pkn^^ (3.1) 



(y) = Y J P n (y)n(n-l)---{n-k + 1) (3.2) 



n=k 



and where P n is the probability of having exactly n particles in the system at time (rapidity) y. The 
hierarchy equation is equivalent to the master equation 

dP 

= -anP n + a(n - l)P n ^ + (3(n + l)(n + 2)P n+2 - (3n(n + l)P n+1 (3.3) 

which requires an elementary 1^2 splitting with transition rate a, a 2 — > vertex having transition 
rate (3 and a 2 —* 1 vertex with transition rate —(3. These later two vertices are certainly not natural in a 
genuine reaction-diffusion process, especially the 2 — > 1 vertex with negative probability. Nevertheless, 
formally, such vertices in a reaction-diffusion context are equivalent to a Reggeon field theory having 
only a triple pomeron coupling. The presence of a 2 — > 1 vertex with a negative transition rate means 
that P n is not guaranteed to be positive but, as we shall see, this does not appear to cause difficulties 
in the situation we shall be considering. 

The presence of a 2 —* vertex means that 

^ = 2PP 2 (3.4) 
dy 

so that there are transitions into the zero particle state. However, there are no transitions out of the 
zero particle state, as is apparent from (|3.3j) . so that ultimately the system will settle down completely 
into the zero particle sector. However, in the situation of interest to us, f3 = a 3 is very small so that 
P 2 at early times is very small and thus the lifetime of the system is very long, of size e 1//2 ° 2 as we 
shall find below. Thus the system will first settle down, rather quickly, into a quasi-stable state and 
then slowly all the probability will appear in the n = sector. This slow transition to the n = state 
will correspond to the tunneling of the 3P Reggeon field theory discovered long ago. 
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To be specific suppose at y = we take 



P n (0) = 5 nl , (3.5) 

that is we start with a single particle. So long as a/ (5 3> 1 probability rapidly flows to higher n states, 
according to ()3.3|) . via 1 — > 2 transitions. When re reaches size a/ (5 the 2 — > and 2 — > 1 transitions 
become important and the system becomes metastable. It is easy to check that 

n-\-l 

*" -\ (i) ' (3-0) 

with (—1)!! = 1, gives = for n > 1 when (j3.6j) is substituted into 1)3.3)1 and furthermore that this 
is the case even when a different value of c is taken for even values of n as compared to odd values 
of re. The n = part of ()3.3|) is given by (|3.4|) which, along with the n = 1 part of 1)3.3)1 . is the only 
departure from ()3.6|) being a fixed point solution. Now starting from ()3.5|) we expect that P n should 
be smooth in n at large values of n. In order to insure this smoothness one must choose the c in ()3.6j) 
to be different for even and odd values of n. A simple calculation shows that if one takes 

c = 1 for odd n, 

2/ ir for even n (3-7) 

then ()3.6j) is smooth in re, for large re, and P n is a normalized probability distribution when a/ (3 ^> 1. 
One can expect the form ()3.6|) to set in at a "time" 

2/0 - 1/2/? (3.8) 

after the start of the evolution of the system from the initial distribution 1)3.5)1 . 

Now as time goes by more and more of the probability of the system will be located in the n = 
state. This means that one must change the normalization of 1)3.6)) to account for the loss of probability 
in the normal re ^ states. This is easily done by including a factor 1 — Po(y) on the right hand side 
of ()3.6)) . One can then, in general, write 1)3.4)1 as 

^-(1-Po) = -2/?P 2 (l-P ) (3.9) 
where P2 is given exactly by 1)3,6)1 . This equation is easily solved as 

1/2 



1 - Po(y) = exp 







V 




V2^ 





e ~ a /2f3 



(3.10) 
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This is exactly the tunneling probability found by Alessandrini et al. jSHl E3 EH1 OH BQI and given in 
(|2.41j) above. 

The factor multiplying — y in the exponent in (j3,10|) is to be identified with otp — 1 with ap the 
pomeron intercept. The high energy scattering amplitude will decrease at large y a result here which 
is easily interpreted as the loss of probability from the active states, n 7^ 0, to the inactive, n = 0, 
state. In the reaction diffusion picture there is no tunneling, only a very slow loss of probability from 
the metastable fixed point configuration of probabilities given by (|3.6|) into the n = state. 



3.2 3P and AP couplings 

The hierarchy equation (see ()2.43j) 



aknW + ak(k - l)n (fc-1) - (3kn {k+1) - (3k(k - l)n (fc) (3.11) 
dy 

along with the corresponding master equation 
dP 

= -anP n + a(n - l)P„_x + [3n(n + l)P n+1 - (3n(n - l)P n (3.12) 

dy 

correspond to a genuine reaction diffusion process with a 1 — ► 2 vertex of strength a and a 2 — > 1 



vertex of strength f3, see e.g. Ref. |22]. In Reggeon field theory language the additional term in (|3.11|) . 
as compared to (|3.1[) . corresponds to a AP coupling. The fixed point solution to (|3.11|) and (|3.12j) is 

p - = ^(i)"w<b ("3) 

and 



Eq. (|3.12j) does not have transitions to the n = state so if one takes a normalized initial condition 
having Pq = then Po(y) will remain equal to zero. When a/f3 3> 1 the fixed point solution is 
very close to a Poisson distribution of probabilities. The fact that there is now a genuine fixed point 
solution to the reaction-diffusion problem guarantees that ap = 1 in the solution to the Reggeon field 
theory, and this indeed is the case as we have seen in Sec. 2.5. 
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3.3 High energy amplitudes 



As discussed earlier in Sec. 2.3 it is not so straightforward to decide how to define a high energy 
scattering amplitude. In Reggeon field theories one has to decide how the pomeron field, the dynamical 
variable in the theory, couples to external particles which are not part of the dynamics of the theory. 
The most common choice, given in (|2.24JI and (|2.25j) . has been to directly couple the pomeron field 
to a source at y = with strength, say, g± and to another source of strength, say, gi at y = Y. In 
the following we shall deal exclusively with the model described in Sec. 3.2 as this seems the more 
interesting case. The 3P case is easily treated by the same methods and the result is as in (|2.41j) . 

Once an s-channel picture of the dynamics is introduced more constraints on the definition of the 
scattering emerge. In the first place the particles which scatter are the same as the particles making 
up the s-channel dynamics. Secondly, the scattering is no longer manifestly boost invariant. For 
example, in the center of mass the scattering looks like a scattering of two highly evolved systems 
while in the rest frame of one of the particles the scattering looks like that of a highly evolved system 
on an elementary quantum. For both (|3.1|) and (|3.11f) boost invariance can be satisfied by taking the 
scattering amplitude to be 

T ( y ) = E 1 a 2k n^ k \Y )n^(Y - Y ) (3.15) 

k=i fc - 

where Yq is arbitrary. Eq. Q3.15|) represents scattering as the sum over k interactions of elementary 
quanta in one, say left-moving, system evolved to rapidity Y — Yq and another, right-moving, system 
evolved to rapidity Yq. The elementary interactions are of strength c? and T is independent of Yq if 
a 3 = p. The use of factorial moments seems necessary in (|3.15|) in order to get boost invariance. At 
large y values rS k \y) approaches the fixed point solution given in (|3.14fl . independent of the initial 
condition, in the case of a genuine reaction diffusion process corresponding to 3P and 4P interactions. 
T is most easily evaluated by taking Yq = in which case 

T(Y) = a 2 n (1 \Y) (3.16) 

which has the large Y limit 

T (Y) F ~ a 3 /P - 1 + e~^ 2 (3.17) 

where we have taken a 3 = j3 in the right-most expression in (|3.17j) . This agrees with the result found 
from Reggeon field theory in H2.56j) . The fact that T is very near one at large Y and small a follows 
from the fixed point properties of the theory along with boost invariance and is independent of the 
initial condition for the system so long as the initial states are normalized to unit probability. 
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Now let us evaluate the scattering amplitude in the eikonal model, repeating the derivation given 
in Sec. 2.5 but from a different perspective. It is still convenient to use reaction-diffusion dynamics, 
but we shall define T in the way that has been traditional in Reggeon field theories. That is, we take, 
corresponding to Q2.25)l . 



T (Y) = £ ( 5 2«) fc n«(y) (3.18) 

fc=l 

where satisfies the hierarchy equations 1)3. 11|) of the Reggeon field theory. The initial condition 
for is 

n^(0) = ( gi /a) k (3.19) 

expressing the eikonal model of pomeron interactions with the external particle (source) at y = 0. 
Now ()3.19)1 requires 

Pn(0) = -.{ gi /a) n e-^ a . (3.20) 
n! 

We note that 

oo 

£p n (0) = l- e "*/« (3.21) 

n=l 

so that the initial state is not properly normalized from the s-channel (reaction-diffusion) point of 
view. Eq. ()3.21)1 gives 

n (k) (Y) y ~ (l/a 2 f (1 - e-«/«) ^ _ ^ (3.22) 
which, when used in 1)3.18)1 . gives, as in 1)2.55)1 . 

T(Y) y ~ (1 - e-^ a )(l - e-»/«) - _ ^ (3.23) 

or 

T(X) yZToo (1 " e" 52/a )(l - e^/ Q ) (3.24) 
where we have dropped small terms, of size e -1 /" 2 , in obtaining 1)3.24)) which is the same as given in 

Refs. |2g||22|[X2|. 

Eq. 1)3.24)1 expresses "gray" rather than "black" scattering. From the Reggeon field theory point 
of view this result is rather mysterious; why should a strongly interacting scattering not saturate 
unitarity at high energy? From the reaction-diffusion (s-channel) picture the factors in ()3.24j) just 
represent the normalization of the states (see 1)3.21)1 ) involved in the scattering. 
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3.4 Generating functions 



For completeness, let us finish with a brief discussion of a generating function, Z(y, v), for occupation 
number probabilities, P n (y), in the reaction diffusion process, 



Z(y,v) = "£P n (y)v n , 



n=l 



where v is an auxiliary variable. The probabilities of n-particle states may be expressed as 



and the factorial moments read 



1_ d»Z{y,v) 



(fe) , _ d*Z{y,v) 



v=0 



« w (y) 



Q v k 



(3.25) 



(3.26) 



(3.27) 



v=l 



The generating function is intimately connected to the wave function of the corresponding RFT, given 
by (j220J). Using (|2~TT|) . (|2~20|) . the identification (u k ) y = n^ k \y) and the relation in (|3~27|) we obtain 

1 d k Z(y,v) 



fc=i 



k\ dv k 



v=l 



which may be simplified to 



9(y,q) = Z(y,l)-Z(y,l-aq). 
It is straightforward to verify that the inverse relation reads 

Z(y,v) = V(y,l/u) -V(y,(l-v)/a). 



(3.28) 



(3.29) 



(3.30) 
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